BUTP-95/18 

Fixed point action and topology in the CP model 

Rudolf Burkhalter 
Institut fiir theoretische Physik, Universitdt Bern, 
Sidlerstrasse 5, CH-3012 Bern, Switzerland. 
(December 1995) 

Abstract 

We define a fixed point action in two-dimensional lattice CP^~^ models. 

The fixed point action is a classical perfect lattice action, which is expected 

to show strongly reduced cutoff effects in numerical simulations. Furthermore, 

the action has scale-invariant instanton solutions, which enables us to define a 

correct topological charge without topological defects. Using a parametriza- 

tion of the fixed point action for the CP^ model in a Monte Carlo simulation, 

we study the topological susceptibility. 
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I. INTRODUCTION 



One possible way of regularizing a continuum quantum field theory is to introduce a 
lattice as a UV regulator. An additional benefit of this method is that it opens the door 
to computer simulations. However, by naively discretizing physical observables, the correct 
values are only obtained in the continuum limit, when the lattice spacing is going to zero. At 
finite lattice spacing the lattice induces systematic errors (cutoff effects). In order to remove 
these cutoff effects one has to introduce finer and finer lattices, and finally extrapolate the 
calculated quantities to the continuum limit. This delicate limit is a major difficulty in 
extracting continuum physics from the lattice. 

It is possible to circumvent this problem by using Wilson's renormalization group theory. 
If one constructs the lattice action and all operators in a form that corresponds to the 
renormalized trajectory, one obtains results that do not depend on the lattice spacing. It 
has been shown in a series of papers recently published that it is possible to construct 
a fixed point action for asymptotically free theories. This action is the fixed point (FP) of 
an exact renormalization group (RG) transformation and as such may be taken as a first 
approximation to the renormalized trajectory. 

As a pilot project a local parametrization of the FP action was constructed in the 0(3) 
nonlinear a model, and used in numerical simulations [|I|]. The result was very promising: 
Although the FP action is perfect strictly only classically, no cutoff effects were seen even 
at small correlation lengths (^ ~ 3). 

In a subsequent paper ^ a FP topological charge was proposed. It was shown, that the 
combination FP action and FP charge has no topological defects. The FP action has the 
correct value for instanton solutions, and does not depend on the scale of the instanton. 
Hence it admits stable instanton solutions on the lattice. This is in contrast to the standard 
lattice action, whose value depends on the instanton scale, and which suffers from disloca- 
tions 0]. These dislocations were suspected to be responsible for the non scaling behavior of 
the topological susceptibility in the 0(3) nonlinear a model. However, using parametriza- 
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tions of the FP action and the FP topological charge in a Monte Carlo simulation, even with 
the absence of dislocations a strong violation of scaling of the topological susceptibility was 
established 0. This indicates that the topological susceptibility is not a physical quantity 
in the 0(3) nonlinear a model. This observation is supported by results of semiclassical 
approximations where the instanton size distribution is divergent for small instanton 
sizes. 

In CP^~^ models with N > 2, however, semiclassical approximation indicates that there 
is no dominance of small instantons. Furthermore, for > 3 dislocations are suppressed 
even for the standard lattice action It should, therefore, be possible, to determine 
unambiguously the topological susceptibility in these models. On the other hand, there 
were several recent determinations of the topological susceptibility in the CP'^ model, using 
different discretizations of the action and the topological charge [[7|-p!T[]; but the results of 
these determinations are partly in plain contradiction with each other. The situation is 
therefore by no means clear. We suggest to use the same concepts and methods in order to 
determine the topological susceptibility as was proposed in Ref. 0]. 

This paper is organized as follows. In Sec. |ITT| we construct a classical perfect lattice 
action for the CP^ model along the same lines as in Ref. We define a RG transformation 
and determine the FP with analytical and numerical methods. We give a parametrization 
of the FP action, that works reasonably well even for coarse grained fields. In Sec. |V| we 
define a FP topological charge along the lines of in Ref. |^. In order to use the FP charge 
in numerical simulations, we construct a parametrization of the dependence of the fine field 
on a coarse input field. In Sec. we discuss the topological susceptibility and the influence 
of a lattice regularization on its measurement. Finally, we present the results of numerical 
simulations using the FP action and the FP charge. 
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II. MODELS 



Two-Dimensional CP^~^ models are an important testing ground for methods in quan- 
tum field theories because of their similarities with four-dimensional non-Abelian gauge 
theories. Important common properties are asymptotic freedom, dynamical mass genera- 
tion, confinement of non-gauge-invariant states and a nontrivial topology. Despite these 
common features, two-dimensional spin models are much easier to handle both analytically 
and numerically. 

The CP^^^ model (in the continuum) consists of a A^-component complex spin field 
z^{x) which satisfies the constraint z{x) ■ z{x)=l. The action has a global SU(A^) symmetry 
and a local U{1) (gauge) symmetry. The gauge field, however, is not independent and 
can be expressed in the basic field z{x). Only after quantization, the gauge field becomes 
dynamical, and gives rise to a confining potential, as is explicitly seen in the large- limit 
|T^. In the following it will sometimes prove useful, not to use the basic field z{x), but the 



gauge invariant composite operator P(x) = z(x) (g> z(x). For N = 2 the CP^^^ model is 
equivalent to the 0(3) nonlinear a model and P{x) can be written in terms of the usual 
0(3) covariant spins S(a;) as P{x) = i(l] + a ■ S{x)). Here a are the usual Pauli matrices. 
In the notation with the composite operators the continuum action for CP^^^ assumes the 
form 

iV/?Aont = ^ / d^^tl {d^P{x)d^P{x)} . (1) 

One has several possibilities to discretize the continuum theory, and to put it on a lattice 
|T3| . If one approaches the continuum limit, each of it will give the same results. We choose 
for this study a discretization, which does not use an explicit gauge field. The standard 
lattice action without gauge fields assumes the form 

N(3As^ = AT/?^ {l - \Zn+i,Zn\^} . (2) 

For N = 2 this goes over to the standard lattice action of the 0(3) nonlinear a model, 
which is not the case if one would use a formulation with an explicit gauge field. The 



parametrization of the FP action will be a generalization of the standard action, taking into 
account couplings between two spins, that are more distant than nearest neighbor and also 
multispin couplings between up to four spins. 



III. FIXED POINT ACTION 

A. Equation for the fixed point action 

We consider the CP^^^ model on a two-dimensional square lattice with variables Zn at 
each lattice site n. Then we perform an exact RG transformation. For this we split the 
lattice into 2x2 blocks hb and with each block we associate a block spin (ng- We define 
the (gauge invariant) RG transformation as an averaging over the fine spins in one block, 
that has the form 

^-Np'A'iO = f ^(^) exp + nN/3j2T. ICn^zA . (3) 



Here, k is a free parameter of the RG transformation, and the normalizing factor 7V(z) 
assures, that the partition function does not change under the RG transformation 

M{zy' = j exp Ln^Y. E \CngZn?] ■ (4) 

The form of the RG transformation in Eq. (^) was chosen, such that for the CP^ model it 
corresponds to the one used for the 0(3) nonlinear a model in Ref. [|1|. In order to compute 
the normalizing factor N'{z), consider the unitary matrix 

= Y Zn® Zn- (5) 

It is just the sum of the composite operators P„ at every site of a block. Let the Hermitian 
matrix M„g have eigenvalues A*„^ and eigenvectors w^^. Because M^g is gauge invariant, 
also its eigenvalues and eigenvectors do not depend on the gauge. Using these eigenvalues 
and eigenvectors, we can rewrite one term in Eq. (^ 

E Kn.^nP = CngM^BCng = KbKubKJ ■ (6) 
news i 



In the limit (3 oo we perform a saddle point approximation. We get to lowest order 
Af{z)^^ = Gxp{KNPJ2nB -^ns)' wliere Xnj^ is the largest eigenvalue of M„g. The FP of the 
RG transformation is determined in this limit by the implicit equation (FP equation) 

Ap(C) = mill { Ap(^) + T{(,z)}, (7) 
with the transformation kernel 

T{C,z) = ^Y.(^^n,- E \Cn,-Zn\A. (8) 
ub \ n^UB / 

The kernel T((^, z) is gauge invariant in a strong sense: It does not change under independent 
gauge transformations on either the coarse or on the fine spins. This is more than would 
be needed: in a general case it would be sufficient, that the kernel is gauge invariant under 
a combined gauge transformation on the coarse and the fine spins. The gauge invariance of 
the kernel assures, that if the action ^pp(2;) of the fine field is gauge invariant, then also the 
action ^pp(C) of the coarse field after the RG transformation is gauge invariant. 

The FP equation (|^) determines the value of the FP action for a given input configuration 
{C}- One may solve this equation iteratively, leading to a minimization on a multigrid of 
lattice configurations with the configuration {C} on the coarsest level, and with k finer 
configurations {z^^^} on successive levels: 

^(^•)(C)= min U(°)(2(^)) + T(z(^-i),2('=)) + ... + r(C,2('))|. (9) 

On each successive level the spin configurations become smoother and smoother, hence one 
may choose for the action A^'^\z^'^^) on the finest configuration {z^'^^} any lattice discretiza- 
tion of the continuum action. The FP action ^fp(C) is then obtained in the limit k ^ oo 
of A^''^{C). For practical purposes, however, only a few levels are needed, and, starting from 
the standard action on the lowest level, the FP value is reached soon. This iterative method 
can be used to solve the FP equation numerically. One can, however, make some important 
statements even without explicitly solving Eq. (|^). This will be done in the next section. 
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B. Classical Solutions 



As in the 0(3) nonlinear a model we can show the following important statement con- 
cerning classical solutions (e.g. instanton solutions). 

Statement: If the coarse configuration {(^} satisfies the FP classical equations of motion 
(i.e. the classical equations corresponding to ^fp), and therefore is a local minimum of 
v4.Fp(C), then the configuration {-2(C)} on the fine lattice, which minimizes the right hand- 
side of Eq. (|^, satisfies the equations of motion as well. In addition, the value of the action 
remains unchanged: Ay-p{.z{Q) = ^fp(C)- 

Proof: Since {(} is a solution of the FP classical equations of motion, it is a stationary 
point of ^Fp(C)) it satisfies 6A-e-p/6( = 0. The transformation kernel is positive: 

Xrs = - CnsMnsCns > 0- (lO) 

The minimizing configuration {z} must fulfill 7^^ = 0, otherwise we could find another 
coarse configuration {(}, that fulfills this equation and lowers therefore the value of ^Fp(C)j 
which would be in contradiction to the assumption, that {(} is a stationary point. Hence 
we have a (albeit not unique) relation between the given coarse field {(} and the minimizing 
fine field {z{C)} 

MubCub = ^ubCub-i (11) 

in other words, is the eigenvector corresponding to the largest eigenvalue of the matrix 
MuB defined in Eq. (|^). T(C, 2;) is zero (that means, at its absolute minimum) for the 
minimizing fine field {z{C)}, hence {2;(C)} is also a stationary point of ^fp and {-2(C)} 
satisfies the classical FP equations of motion. Furthermore, the value of the action is the 

same: ^fp(-2(C)) = '^fp(C) q-e.d. 

Note, that the reverse statement is not always true. If the fine configuration {z} is a 
solution of the equations of motion, then the coarse configuration obtained from Eq. (pUD 
is a local minimum, but it need not be the absolute minimum that has to be found at the 
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right-hand side of the FP equation (|^). This mechanism prevents the existence of arbitrarily 
small instantons on the lattice (see below). 

Using the above statement, we can now construct instanton configurations on the lattice 
(see also Ref. ||^ for the corresponding construction in the 0(3) nonlinear a model). We 
consider instanton configurations of the continuum on a torus. On a torus, however, there 



exists no exact one- instanton configuration [|14|. In order to clearly separate cutoff effects 



from finite size effects, we have to fall back upon exact two-instanton configurations. The 
two- instanton configurations of CP^, embedded into CP^~^, have the form 

2 

\ a(w - hi) 



z{w) = M 



T-r <^\w — a 

u+ n —' — —V 



(12) 



w = x + iy, M = (1,0,0,...), t;= (0,1,0, ...), ^ai = ^h 



k 



1=1 1=1 



where o"(?t,) is the Weierstrass a function and the factor J\f ensures the correct normalization. 
The vectors u and v are A^-component CP^~^ vectors which specify the orientation of the 
instantons in "color space". The four complex parameters ai, 02, &i, and 62 specify the size 
of the instantons and their position and orientation on the torus. We choose them in the 
form 

ai = {L/2-p)i, 02 = L/2 + (L/2 + p) z, (13) 
61 = (L/2 + p) i, 62 = L/2 + (L/2 - p) i, 

where p is the instanton size and L the size of the torus. The action ^cont of the con- 
figuration (O) is equal 2 x 27r irrespective of the size p of the instantons. Next we want 



to construct lattice two-instanton configurations of different sizes. We first discretize the 
continuum two-instanton configuration on a very fine lattice, so that the cutoff effects can 
safely be neglected. Then we perform k block transformations using Eq. (pAf ). Under a block 
transformation the size of the instantons is halved. Choosing the number k of block trans- 
formations and the initial size p, one can get any final size p 2"^^ on a coarse lattice. The 
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above statement shows, that the action remains the same, unless the size of the instantons 
is too small, and they fall through the lattice. 

We have numerically performed the above program for the CP^ model, and have mea- 
sured several quantities on the finally blocked configurations. On the coarse configuration 
itself we measured the geometric charge, the standard action, the Symanzik improved action 
and the parametrization of the FP action presented in Sec. [Ill Q . Performing a minimization 



on a multigrid with three finer levels, we measured the exact FP action and on the finest 
level the FP charge (cf. Sec. ITVl) . Because instantons in CP^ are embedded CP^ instan- 
tons, it is clear, that the results presented in Fig. ^ are practically identical to the ones 
obtained in the 0(3) nonlinear a model 0. The results show, that it is possible to obtain a 
parametrization of the FP action that performs very well for instanton configurations down 
to the smallest possible size on the lattice. In contrast, the standard action and also the 
Symanzik improved action perform quite bad, especially for small instantons. It is worth 
mentioning, that the range of instanton sizes, where the geometric charge differs from the 
FP charge, is quite narrow. For instanton configurations, the geometric charge seems to 
be as good as the FP charge. For general configurations created in a Monte Carlo (MC) 
simulation, however, there is a noticeable difference between the two charges (cf. Sec. |VI| ). 

C. Parametrization of the FP action 

If we want to use the FP action in practical numerical simulations, we must construct a 
parametrization, that consists of a not too large set of local operators. Our parametrization 
has the form 

'4pp''(^) = -^IIp('^)^n,n+r + H coupliug X products of , (14) 

n,r ni,nj,... 

where 9rn,nj = arccos Zn^Zn^ ) is the angle between two spins. Note, that in CP^~^ models 
the maximal angle between two spins is 7r/2; z and —z is the same spin, different only 
by a gauge transformation. There are two reasons why it is useful to use the angle 6'„._„j. 

9 



instead of (^1 — p^: For solutions of the equations of motion of the form z{n) = 

{cos6no, sin^^riO) 0,...), the 6 dependence of the action is exactly 6"^. Moreover, if one 
rotates a single spin in a trivial background with all spins pointing in the same direction, 
the 6 dependence of the action is more like 6"^ rather than 1 — cos^(^). Figure |^ shows how 
well the minimized FP action is approximated even by the two first terms (nearest neighbor, 
diagonal) of the lowest order of the parametrization in Eq. (|T^). In contrast, the standard 
action performs very bad for large angles. 

We can calculate the coefficients p(r) of the lowest order analytically. The coefficients of 
the higher orders can be determined in a numerical fitting procedure. 

1. Lowest order, determination of p 

The analytic result for the lowest order term in the FP action will be valid for all A^, 
as long as > 1. Consider on the coarse level a smooth configuration, which is weakly 
fiuctuating around the first direction 



where has (A^ — 1) components, and ^ 1. With this choice we have fixed the 

gauge by imposing, that the first component is real and positive. Then the minimizing fine 
configuration will also fluctuate around the first direction and, using the same gauge fixing 
prescription, we make the ansatz 

/. /l _ IT/ |2\ 

(16) 



l-IK 



with \Vn\ -C 1. Putting these expansions and the parametrization of the action (14) in the 
FP equation (j^) and keeping only terms up to quadratic order in X and V, we get the 
equation 



^}^\IY1 p{r)Re{VnVn+r) + 2k ^ 

i^J I ^ n,r ub 



10 



riGnB 



2 



(17) 



This equation can be solved for p most easily by using the same technique as described in 
Ref. []I|. Taking into account the fact, that for the CP^ model the RG transformation goes 
over to the one used for the 0(3) nonlinear a model, it is not astonishing, that the resulting 
p is exactly the same as in Ref. [0]. Consequently, in order to have a most local action, we 
will also choose the free parameter = 2 in the RG transformation Eq. (j^). 



2. Higher orders for the CP^ model 

In order to have a reasonably good parametrization of the lattice action which performs 
well also for coarse configurations, we include higher order terms in the parametrization (|T4|). 
The coefficients of these higher order terms will depend on the chosen model, and are 
different for each N . We determined them numerically for the CP^ model. To do this, 
we first produced about 300 configurations of lattice size 3 with a Monte Carlo program 
using the standard action. The configurations ranged from ones with small actions and 
fiuctuations with small amplitudes to strongly fiuctuating ones. We also included some two- 
instanton configurations with radii of the order of one lattice spacing (cf. Sec. [Ill B] ). For 
every configuration we calculated the value of the FP action by minimizing the FP equation 
on a multigrid. Then we determined the coefficients of 30 higher order operators chosen 
because of their locality. Their value was determined so that the average difference between 
minimized and parametrized FP action was minimal. 

Figure ^ illustrates the quality of the parametrization. The relative deviation from the 
minimized FP action is only large for very coarse configurations with large actions. For 
comparison, we also plot the standard action for the same configurations. Note, that here 
the relative deviation of the standard action from the FP action is not small, even for 
configurations with small amplitude fiuctuations. 



The resulting 32 coefficients of the parametrization ( [I^ ) are given with a graphical no- 
tation of the corresponding operators in Table |. A line with two dots • — • means that 
the angle 6"^^ between the two spins at positions rii and 77,2 of the dots enters into the 
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parametrization. A graph consisting of several lines represents just the multiplication of the 
corresponding angles. Coefficients No. 1 and No. 4 are the only analytically calculated ones, 
all the other couplings were determined by the fitting procedure. The locality of the action is 
expressed in the fact, that it is possible to construct a good parametrization with operators 
consisting of spins, that lie only within a 2 x 2 section of the whole lattice. On the other 
hand, the parametrization seems to be not as local as the one found in the 0(3) nonlinear a 
model. Some couplings of operators of higher order have quite a large value. Nevertheless, 
one should not give too much attention to the actual value of the couplings. A large coupling 
for a higher order operator does not mean, that the action is less local. Some combinations 
of couplings can be changed without much affecting the quality of the parametrization. In 
principle one would not need to keep all 32 operators in the parametrization, we just did 
so for reasons of completeness. Using 32 operators for an action in a MC simulation may 
seem to be a lot. However, compared with the expected benefit of using such an action, the 
additional computational effort needed seems to be reasonable. 

IV. FIXED POINT TOPOLOGICAL CHARGE 

A feature of CP^~^ models common with gauge field models is the existence of topologi- 
cally nontrivial solutions. In the continuum theory the topological charge Q may be defined 
as the integral 

Q = ^ 1 rf'a; tr \P{x)d^P{x)d,P{x)\ . (18) 

The action is related to the charge through the inequahty 

Aont>27r|g|. (19) 

For instantons the equality holds. They minimize the action, and are therefore solutions of 
the equations of motion. 

On the lattice, however, this concept breaks down because continuity is lost. In the 
continuum the topological sectors are clearly separated, but on the lattice one may continu- 
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ously transform a field from one topological sector to another. Furthermore, when using the 
standard lattice action, the scale invariance of instanton solutions is violated. The action 
decreases with decreasing instanton size, and configurations that violate the inequality (|19|) 
- so called "dislocations" - are possible. 

In Refs. it was shown that the FP action admits stable instanton solutions. Fur- 
thermore, it was shown how to construct a correct charge - the FP charge - which does not 
admit dislocations. In this section we will proceed along the same lines as in these works, 
in order to construct the FP charge for the CP^ model. 



A. Definition of tlie FP ciiarge 

We define the FP charge by means of the iterated solution of the FP equation (^) . Under 
a RG transformation in the limit (3 oo, an operator 0{z) transforms into 0'{Q on the 
coarse lattice as 

O\0 = 0{z{C)), (20) 

where the spin configuration {-2(C)} is the solution of the FP equation (|^). The FP topo- 
logical charge is obtained in the limit of infinitely many RG transformations: 

Qpp(C)= limQ(^('^)(C)), (21) 

fe— >oo 

where {z^^^*} is the solution of the iterated FP equation (P) on the lowest level in a level 
multigrid. As the configurations get smoother on each successive level, one may choose for 
the charge Q on the lowest level any lattice discretization of the topological charge. In this 
paper we will use the geometric charge |T^, because it is stable against small variations of 



the field, if the field is sufficiently smooth, and because it always gives an integer number. 

In Ref. [0 it was shown, that the combination FP action and FP topological charge 
always obeys the inequality 

Ap(C)>27r|Qpp(C)|. (22) 
13 



Hence there are no dislocations. 

In numerical simulations it is very time consuming to minimize a multigrid for every 
configuration. One needs a parametrization of the solution {2;(C)} of the FP equation. This 
will be done in the next section. 



B. Fixed point field 

The FP field is the fine field z^''^ in the multigrid solution of the iterated FP equation @ 
as k goes to infinity. If the functional dependence of the solution on the first fine level 2;^^^ 
on ( is known, the FP field can be evaluated by iteration. Below we give an expression for 
the field z^-^^ = z^^\(). (In the limit k 00 the solution {-s;*-^^} of the iterated FP equation 
is identical to the solution {z} of the FP equation.) 

We make the expansions (0) and (|I6]) for the coarse and the fine configuration. Then 
the solution of the FP equation leads in lowest order to the relation 

Vn = Yl "('^' ^B)Xng (23) 
with a the same as for the 0(3) cr model [1T1-0|,|16 



To obtain this result, we have fixed the gauge in a specific way. However, we want 
a relation between {z} and {C}, that does not depend on the choice of the gauge fixing 
prescription. Furthermore, for coarse configurations we have to include higher order terms, 
that parametrize the dependence of {z} on {(}. For this we construct the gauge invariant 
composite operators Png = Cub ® Cns on the coarse lattice. On the fine lattice we build 
matrices Qn by, summing over the coarse composite operators including the next to leading 
order terms: 

Qn = Y.a{n,nB)PnB+ ^(^' '^b, m^, "^b)^^,^^^™^' (24) 

n_B "B 

and define the fine field variable Zn as the eigenvector of Qn with largest eigenvalue. In next 
order enters the angle 9'^^ between the coarse spins at sites niB and m'^, respectively. 
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In order to determine the coefficients (3, we used the same ~ 300 configurations hke for the 
parametrization of the action, minimized the FP equation (|^ and stored the resulting fine 
lattices. The coefficients were then determined by minimizing the difference between the 
minimized fine spins and the parametrization (0). 

The numerical values of the coefficients a and (3 are given in Table ^ with a graphical 
notation of the corresponding operators. We chose a set of 17 operators mainly because of 
their compactness. Numbers 1-6 are the analytically determined coefficients a, Nos. 7-17 
are the numerically determined coefficients (3. The meaning of the graphical notation of 
the operators is the following: The dashed lines represent a 3 x 3 section of the coarse 
lattice grid. The cross + in between indicates the position n of the matrix Qn in Eq. (^4|). 
The little square □ denotes the position of the coarse composite operator P„g , the two 
connected dots • — • are the positions mg and of the spins whose angle O"^^ ,^/^ enters into 
the parametrization. Graphs obtained by trivial symmetry transformations are not drawn 
separately. 



V. TOPOLOGICAL SUSCEPTIBILITY 

The topological susceptibility is defined as the ratio 

X. ^ (25) 

where Q is the topological charge and V is the space-time volume. In CP^^^ models it is 
a dimension two quantity, that vanishes to all orders in the weak coupling expansion. From 
the perturbative renormalization group it is expected to scale according to the two loop (3 
function 

/I 

Xt oc (2^^^) exp(-7riV/3) {f3 ^ 00). (26) 

In order to check the continuum limit, a scaling behavior is more important to observe 
than the above asymptotic scaling. One additionally measures a second quantity, e.g., 
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the correlation length ^, and builds the dimensionless product Xt^^; which should go to a 
constant in the continuum limit ,^^00. 

In general, cutoff effects can originate from two sources: from the discretization of the 
action and from the discretization of an operator. The topological charge is an operator, 
that strongly exhibits lattice cutoff effects, as was demonstrated in Sec. p.llB|: instantons 



with a radius smaller than about 0.7a fall through the lattice and are lost in a Monte Carlo 
simulation. This fraction is large at correlation lengths of the same order as the lattice 
spacing, and gets smaller as the correlation length grows. One therefore expects a scaling 
violation, which is determined by the small instanton size distribution. 

One may estimate the behavior of the expected scaling violation by using the results of 
a semiclassical expansion, and performing a kind of dilute instanton gas approximation. In 
CP'^ the probability density to find a field configuration with topological charge Q = 1 in a 
sphere with radius i? is 

D^{R) = K,{4) Ai,R' ^^(A^gi?) [R ^ 0) (27) 

where A^s is the perturbative A parameter and Ki{4) = 36.2995 results from the integration 
over the instanton parameters. This result is valid for small volumes, so we can use it to 
estimate the probability of loosing a charge-one configuration that falls through a small 
lattice mesh. 

We denote the charge measured on the lattice with Q and the lost charge with q. If we 
assume that Q and q are independent, then the topological susceptibility measured on the 
lattice is 

(Q') ..cent 

Xt = - Xt — 

We make a kind of dilute gas approximation by assuming that within each lattice mesh there 
can be independently an instanton or an antiinstanton. The result is 

^ = 2Di. (29) 
16 



Identifying R with the minimal size p = ca oi an instanton on the lattice (c ^ 0.7) provides 
^lat ^ ^cont _ 2;^^(4) a4 ^c^a^ \n\A^sCa), (a 0) (30) 

which can be compared with the results of a Monte Carlo simulation. 

VI. NUMERICAL RESULTS 

We performed Monte Carlo simulations in the CP^ model using the parametrized FP 
action given in Eq. (p!^ . In order to reduce critical slowing down, we implemented a hybrid 
overrelaxation algorithm similar to the one described in Ref. 0. The mass gap m = 
was obtained from the long distance fall off of the correlation function projected to zero 
spatial momentum. We made a minimal fit to the correlation function with the function 
c{e~"^^ _l_ g-m{L-a:)-j -^^ ^.j^g interval x G [xmin, L/2] for different Xmin- The value of m is taken 
at an expected plateau at Xmm ~ ^• 

In order to avoid finite size effects, one has to make measurements in large enough 
volumes. We made measurements for all f3 values in volumes with the ratio L/C, — 5.5 — 6 
kept approximately constant. These volumes are usually large enough to totally avoid 
any finite size effects. In CP^^^ models, however, the basic z particles are subject to a 
confining potential. Thus the resulting bound states may have a radius that is larger than 
the correlation length. For this reason we also made measurements for some (3 values in 
even larger volumes (with ratios L/C, up to 12). 

The masses that are determined in the "small" volume L/C — 5.6 are about 3% below 
the ones measured in "infinite" volume, as can be seen in Table pT[ The magnitude of this 



finite size effect fits nicely with the one observed in Refs. P,p!7|. Furthermore, we observed, 
that in the large volumes the determination of the mass gap showed a nice plateau behavior. 
Such a plateau was sometimes not very clearly seen in the small volumes. Nevertheless, we 
may use the results obtained in the small volume in order to look for a scaling behavior in 
XtC"^- The actual values are spoiled by finite-size effects, but since L/C is constant, these 
effects are the same for every point. 
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We checked asymptotic scaling of the mass gap according to the perturbative lattice scale 
in two-loop approximation 

= [-^rcNpj exp [-^rrNpj . (31) 

In Fig. ^ we show the ratio mj bSj^ versus the correlation length, with masses determined 
in "small" volumes Lj^ 5.6 and also some measurements in "large" volumes. The most 
striking observation is that the ratio approaches a constant value m/ A.\ = 8.1(1) (this value 
is obtained in "large" volumes), and that its final value is attained already at a quite short 
correlation length. Such a precocious asymptotic scaling has not been stated for other lattice 
actions. With the standard action, for example, one does not even see an asymptotic scaling 
at correlation lengths ^ ~ 50 PJTO|. Note that the asymptotic scaling behavior for the FP 
action is a purely phenomenological observation that has not been expected theoretically. 
The FP action was constructed nonperturbatively in order to reduce cutoff effects. That the 
mass gap scales with the perturbative A parameter is just an additional, unexpected benefit. 
We furthermore observe, that the value of the ratio mj is remarkably small. From this 
we conclude, that App is much closer to the continuum scale than Ast- 

We determined the topological charge, using both the geometric definition and the defi- 



nition of the FP charge given in Sec. [I^. For the measurement of the FP charge we used the 
geometric charge on a finer lattice of the multigrid with the Monte Carlo generated lattice 
as coarsest level. In order to determine the configuration z{C^ on the first finer level, one can 
either minimize the FP equation (which is very time consuming) or use the parametrization 
of the dependence on ^, given in Sec. |IV B| . We denote the corresponding charges Qcoarse 



for the geometric charge and, e.g., Qpari. level for the charge measured on the first finer level 
using the parametrization of the fine field. 

For some /3 values we compared the results of using the parametrization on a finer level 
and of minimizing on a multigrid. The results were found to be consistent within the 
statistical errors, as is reported in Table |T|. This shows that the parametrization performs 
well for typical configurations occurring in a Monte Carlo simulation. We wanted to be sure. 



that it is sufficient to measure the charge only on the ffist finer leveL As a test we calculated 
for one /3 value also the charge on the second finer level. The result presented in Table |T| 
shows, that the values on the first and the second finer level were found to be consistent 
within the statistical errors. This means that the process of going to a lower level is already 
stable at that stage, and that it is sufficient to calculate the fine field only on the first finer 
level. This is not unexpected, since we use a geometric definition of the topological charge 
at a lower level. Therefore, we always get an integer number for the charge, and the charge 
has to change abruptly when going to a finer level. The field on the first finer level is then 
(in most cases) smooth enough to yield already the FP topological charge. 

Table |T| shows the effect of using different definitions of the lattice topological charge. 
On the coarse level the value of (Qcoarse) higher than the one obtained with the FP charge 
(on the first finer level). This can be explained with dislocations which contribute to the 
geometric charge. They have the effect that one overestimates the amount of topological 
excitation. 

In Fig. ^ we show the results of the scaling test for the dimensionless quantity Xf^^; 
measured in volumes with a ratio L/C, — 5.6. One clearly sees the expected raise at small 
correlation lengths, which is due to lost small instantons in this region. At correlation length 
^ ~ 10 this effect is already saturated. The measurements in "small" volumes using the FP 
charge (i.e., Qpari. level) show the expected scaling plateau at a value of about Xt^^ — 0.071 
(c.f. Fig. The measurements using the geometric definition of the topological charge 
give somewhat larger values but they seem to converge to the value obtained from the FP 
topological charge for larger correlation lengths - these measurements do not show the same 
nice scaling behavior as the measurements using the FP charge. This is also what is to 
be expected in the CP^ model if dislocations contribute to the geometric charge. In the 
continuum limit the effect of these dislocations is suppressed, but they still contribute at 
the considered correlation lengths. 

We can now use the four measurements in "large" volumes, with Qpar i. levei for the topo- 
logical charge and perform a fit to the data with Eq. (|30D , in order to extrapolate to the 
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continuum value. We obtain the value Xt^^ = 0.070(2). 

Let us compare our numerical value with previous determinations and with results from 
the large expansion |Q. To leading order in one gets 

Xte = ^ + 0(l/iV5/3) ~ 0.06 (for iV = 4). (32) 

This value is comparable with the value we got at the scaling plateau. However, the cor- 
rection to the leading order is large in Eq. (^) and the agreement with numerical results 
occurs at chance at = 4. The numerical results of Ref. [§ (who use actions with an 
explicit gauge field) - quoting Xt^"^ — 0.06 with an uncertainty of 10 — 20% - are in quite 
good correspondence with our result. Measurements using the standard action and the ge- 
ometric definition of the topological charge (Refs. P,p!0[]) lead to a value roughly twice as 
large as ours, but with large scaling violations at large correlation lengths - this is probably 
the effect of dislocations which still contribute at the considered correlation lengths. 



VII. CONCLUSIONS 

In CP^~^ models it is possible to define a FP action and a FP topological charge. The 
FP action is the fixed point of an exact RG transformation. It is a classically perfect action 
and possesses scale invariant instanton solutions. The definition of the FP topological charge 
is based on the FP field operator. Both, the FP action and the FP field can be evaluated 
to any precision desired on a sufficiently large multigrid. It can be shown that the FP 
topological charge together with the FP action have no lattice defects. 

It is profitable to use the FP action and the FP topological charge in numerical simu- 
lations. For this purpose, we have parametrized for the CP^ model the FP action and the 
field solution of the FP equation which is iterated to obtain the FP field. We used these 
parametrizations in MC simulations of the CP^ model. We find two main results. First, 
the mass gap unexpectedly scales according to the perturbative lattice scale. Second, the 
dimensionless quantity Xt^"^ raises as expected at small correlation lengths, and reaches a 
scaling plateau at Xt^^ = 0.070(2). 
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For the future, it would be interesting and profitable to use the parametrized FP action 
for further MC simulations, for example in order to investigate the spectrum of CP^~^ 
models. 

ACKNOWLEDGMENTS 

I would hke to thank P. Hasenfratz, F. Niedermayer, M. Blatter, U.-J. Wiese, A. Papa, and 
P. Kunszt for helpful discussions. This work was supported in part by the Schweizerischer 
Nationalfonds. 



21 



REFERENCES 



[1] P. Hasenfratz and F. Niedermayer, Nucl. Phys. B414, 785 (1994), |hep-lat/9308Q0^ ; 



P. Hasenfratz, Nucl. Phys. B (Proc. Suppl.) 34B, 3 (1994); F. Niedermayer, ibtd., 513. 

[2] T. DeGrand, A. Hasenfratz, P. Hasenfratz, F. Niedermayer, Nucl. Phys. B454, 587 
(1995), |hep-lat/9506030| ; B454, 615 (1995), hep-lat/950603l| ; F. Farchioni, P. Hasen- 



fratz, F. Niedermayer, A. Papa, ibid., 638, |hep-lat / 9506032| . 

[3] M. Blatter, R. Burkhalter, P. Hasenfratz and F. Niedermayer, Nucl. Phys. B (Proc 
Suppl.) 42B 799 (1995), |hep-lat/9411068| ; Phys. Rev. D 53, 923 (1996), [hep^ 
lat/ 9508028 . 



[4] M. Liischer, Nucl. Phys. B200[FS4], 61 (1982). 

[5] P. Schwab, Phys. Lett. 118B, 373 (1982); 126B, 241 (1983); M. Liischer, Nucl. Phys. 
B205[FS5], 483 (1982). 

[6] A. Jevicki, Nucl. Phys. B127, 125 (1977); D. Fdrster, ibid. B130, 38 (1977); V.A. Fa- 
teev, I.V. Frolov and A.S. Schwarz, ibid. B154, 1 (1979); B. Berg and M. Liischer, 
Comm. Math. Phys. 69, 57 (1979). 

[7] K. Jansen and U.-J. Wiese, Nucl. Phys. B370 762 (1992). 

[8] M. Campostrini, P. Rossi and E. Vicari, Phys. Rev. D 46, 2647 (1992); D 46, 4643 
(1992), |hep-lat/9207032 . 



[9] U. Wolff, Phys. Lett. B 284, 94 (1992), |hep-lat/920500T 



[10] M. Hasenbusch and S. Meyer, Phys. Rev. D 45, 4376 (1992); Phys. Lett. B 299, 293 
(1993). 

[11] A.C. Irving and C. Michael, Nucl. Phys. B371, 521 (1992). 

[12] A. D'Adda, M. Liischer and P. Di Vecchia, Nucl. Phys. B146, 63 (1978). 



22 



[13] P. Di Vecchia, A. Holtkamp, R. Musto, F. Nicodemi and R. Pettorino, Nucl. Phys. 
B190[FS3], 719 (1981). 

[14] J.-L. Richard and A. Rouet, Nucl. Phys. B211, 447 (1983). 

[15] B. Berg and M. Liischer, Nucl. Phys. B190[FS3], 412 (1981). 

[16] T. L. Bell and K. Wilson, Phys. Rev. B 10, 3935 (1974). 

[17] P. Rossi and E. Vicari, Phys. Rev. D 48, 3869 (1992). 



23 



FIGURES 



FIG. 1. Actions and charge of instantons with radii of the order of one lattice spacing. 
FIG. 2. Actions of configurations with one spin rotated against a trivial background. 
FIG. 3. Scatter plot of the actions of the configurations used for the fit. 

FIG. 4. Asymptotic scaling test for the mass in the CP^ model. 

FIG. 5. Scaling test for the topological susceptibility. These are the results measured in "small" 
volumes L/^ ~ 5.6. The larger error bars at some results from the first finer level are due to the 
fact that we performed less measurements where we also measured the FP charge. 
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TABLES 



No. 


Type 


Coupling 


No. 


Type 


Coupling 


No. 


Type 


Coupling 



1 • • -0.61884 2 



-0.05381 3 



0.20023 



-0.19058 5 



-0.01892 6 



-0.06735 



0.01455 8 



-0.25328 9 



0.07099 



10 



0.04334 11 



0.02704 12 



-0.12660 



13 



0.06787 14 



0.18327 15 



A 



0.13297 



16 



-0.28036 17 



• • 

• • 



-0.00174 18 



0.26017 
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0.11006 20 



0.01396 21 



0.09222 



22 



-0.02530 23 



0.52163 24 



-0.05146 



25 



-0.06314 26 



-0.29456 27 



-0.04937 



28 



0.17930 29 



-0.15733 30 



-0.15941 
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31 



0.11217 32 



0.47978 



TABLE I. Couplings used for the FP action. 



No. Type Coeff. No. Type CoefF. No. Type CoefF. 
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17 


li 1 


0.01849 



TABLE II. Coefficients of the parametrization of the fine field. 
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TABLE III. Results of MC simulations. 
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